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In a recent series of papers, Schwinger discussed a process that he caUed the Dynamical Casimir 

Effect. The key essence of this effect is the change in zero-point energy associated with any change 

in a dielectric medium. (In particular, if the change in the dielectric medium is taken to be the 

r-^ ' growth or collapse of a bubble, this effect may have relevance to sonoluminescence.) The kernel 

O^ , of Schwinger's result is that the change in Casimir energy is proportional to the change in volume 

^\ • of the dielectric, plus finite- volume corrections. Other papers have called into question this result, 

claiming that the volume term should actually be discarded, and that the dominant term remaining 

(~| . is proportional to the surface area of the dielectric. In this communication, which is an expansion 

C^ ' of an earlier letter on the same topic, we present a careful and critical review of the relevant 

analyses. We find that the Casimir energy, defined as the change in zero-point energy due to a 

change in the medium, has at leading order a bulk volume dependence. This is in full agreement 

C'i I with Schwinger's result, once the correct physical question is asked. We have nothing new to say 

about sonoluminescence itself. 
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I. SCHWINGER'S CALCULATION 

Several years ago, the late Julian Schwinger wrote a 
series of papers ^^ wherein he calculated the Casimir 
energy released in the collapse of a spherically symmet- 
ric bubble or cavity. (For the evolution of Schwinger's 
ideas on this topic, see y-QlO Using source theory, he 
derived in a simple and elegant way a formula for the 
energy release involved in the collapse. He found from 
a straightforward application of the action principle that 
(for each polarization state) the "dielectric energy, rela- 
tive to the zero energy of the vacuum, [is given] by 
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So the Casimir energy of a uniform dielectric is nega- 
tive". This was then applied by Schwinger as a model 
for sonoluminescence. 

From this, and following Schwinger, we conclude that a 
dielectric slab of material with a spherical vacuum cavity 
of radius R has a higher Casimir energy than the slab of 
material with the vacuum cavity re-filled with material of 
the same dielectric constant. This energy (per polariza- 
tion state) can be computed by introducing an ultraviolet 
momentum cut-off K into the previous expression for the 
energy [| 
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In general, this volume contribution will be the dominant 
contributionJj 

In view of the elegance and simplicity of this result, it 
is natural to ask whether it can also be derived by more 
traditional quantum field theoretic means. Indeed, the 
existence of such a volume contribution is easy to verify 
on general physical grounds: 



^For generic dielectrics, this wave-number cut-ofF is related 
to the high-wave-number asymptotic behavior of the disper- 
sion relation: it is the scale at which the dielectric disper- 
sion relation begins to approach the vacuum dispersion rela- 
tion. For the particular case of sonoluminescence, this wave- 
number cut-off — deduced from the dielectric properties of the 
medium — implies a minimum wavelength for the electromag- 
netic radiation emitted in the collapse of the bubble. 

^The dots denote finite- volume corrections. We shall inves- 
tigate the leading finite-volume term more fully in a separate 
pubfication M. 



(1) For instance, we know that the effective action in 
3-1-1 dimensions contains divergences which range from 
quartic to logarithmic, in addition to finite contributions. 
Furthermore (as is well known), the "cosmological con- 
stant" contribution (the quartic divergence) never van- 
ishes unless the theory has very special symmetries |9(] 
(as for example, in the case of supersymmetric theories). 
Thus energy densities that go as (cut-off)** are generic in 
(3-1-1) dimensions. 

(2) The cut-off dependence of the energy should not 
be considered alarming: Dielectrics are condensed mat- 
ter systems, and as such abound in physically defined 
and physically meaningful cut-offs — everything from the 
plasma frequency to the interatomic spacing may be con- 
sidered as a candidate for the cut-off. The interesting 
physics comes in deciding exactly which cut-off is physi- 
cally relevant in the current situation. 

(3) We can view Schwinger's result in elementary terms 
as simply the difference in zero-point-energies obtained 
by integrating the difference in photon dispersion rela- 
tions over the density of states 
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The dots again denote finite-volume corrections. The 
energy above is multiplied by 2 because unlike eq. (|l|) 
we have already included the contributions from the two 
polarization states. To actually calculate this energy dif- 
ference wc require a suitable physical model for uj{k). 
(4) Alternatively, we could perform an explicit quantum 
field theoretic calculation of the Casimir energy in some 
specific model problem and thereby verify Schwinger's 
result. We have performed such a calculation, and report 
the results in this paper. 

We find full agreement with the calculation and results 
of Schwinger |l| , and disagreement with some subsequent 
papers. For instance, the calculation presented in ref- 
erences [10 -12 1 computes the Casimir energy associated 
with a spherical ball of radius R, dielectric constant ei, 
and permeability /ii, embedded in an infinite medium 
with dielectric constant €2 and permeability ^2- 

That calculation ultimately asserts that the dominant 
contribution is not the volume term but a surface correc- 
tionQwhich is proportional to R^K^{ei — £2)^- We, how- 
ever, have re-analyzed these calculations and find that 



•^ Actually there is some ambiguity even on this point — 
some, but not all, of the comments in those papers could 
be interpreted as suggesting that the surface term should also 
be discarded, and that the dominant term is actually of or- 
der hc/R. This is tantamount to arbitrarily setting to zero 
anything that depends on the cut-off and simply relying on 
naive dimensional analysis to guess the form of the answer. 
For physical dielectrics with physical cut-offs this suggestion 
is not tenable. 



the volume term is in fact present and dominant except 
for very small bubbles. (The present paper is a more 
detailed and complete presentation of a result previously 
announced in |1^.) 

We have re-analyzed this calculation in several differ- 
ent and complementary ways: 

(A) We have extended one method of doing the calcula- 
tion, namely the method of doing a mode sum over the 
difference of zero point energies, to further illuminate 
the physics behind the cut-off by introducing a num- 
ber of simple general and robust arguments. We show 
that generically surface terms do in fact show up as sub- 
dominant corrections to the dominant volume effect. The 
general arguments contained in this section of the paper 
are particularly useful in that we can deal with arbitrary 
shapes and not be limited by requirements of spherical 
symmetry. This material is reported in section ||. 

(B) We have analyzed previous calculations to see where 
they differ from our own. We first clarify the phys- 
ical content of the different configurations whose zero 
point energies are to be subtracted from each other to 
get the Casimir energy. The results can be compactly 
represented as regulated integrals involving Green func- 
tions. Some prior calculations misapplied the subtraction 
scheme required to get rid of the "vacuum" contribution, 
and we can isolate the difference in energy between these 
prior calculations and the correct one. This energy dif- 
ference can be calculated by an elementary application of 
general quantum field theory results and we thus regain 
Schwinger's volume term. This result is again general 
and not limited by requirements of spherical symmetry. 
The qualitative parts of this analysis, including a general 
demonstration of the existence of the volume term, is in 

and numerical coefficients of the volume term 
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are obtained by this shape independent method in sec- 
tion 0. 

(C) For the case of a spherical dielectric ball, one can 
also address the problem by using expansions of Green 
functions in spherical coordinates, following the formal- 
ism of [pi)[-|l^ . We focus first on the volume term of the 
Casimir energy, i.e., on the energy difference between 
some prior calculations and the full calculation. Given 
our previous discussion, if the volume term is consid- 
ered in isolation the use of specific coordinates makes 
the calculation more complicated than it has to be. How- 
ever, we wish to see the same answer emerge using this 
method. The energy difference is given as a sum over a se- 
ries of integrals involving Ricatti-Bessel functions which 
we evaluate exactly using generalized addition theorems 
for spherical Bessel functions. The result, presented in 
section M, is in complete agreement with the quantum 
field theory calculation and our extensions of the mode 
sum calculation, and also in complete agreement with 
Schwinger's original claim ||l|. 

(D) For completeness, we also present the full calculation 
using the general solutions to the electromagnetic fields 
developed in references pO|-p^ using the dyadic formal- 
ism. We compute the energy difference between (Case I) 



an otherwise uniform medium with dielectric constant 62 
and permeability ^2 , but with a spherical cavity of radius 
R containing dielectric constant ei and permeability /ii, 
and (Case II) a completely uniform medium with dielec- 
tric constant £2 and permeability /i2- This energy differ- 
ence is again given as a sum over a series of integrals in- 
volving Ricatti-Bessel functions. Some of these sums can 
be evaluated explicitly while others can only be evaluated 
by using an asymptotic analysis of the type presented 
in |f0|-[l2[ . We verify the existence of both volume and 
sub-dominant surface contributions, and present these re- 
sults in section VI. 



II. MODE SUM OVER ZERO POINT ENERGIES 

We wish here to present a simple derivation of 
Schwinger's result and to make some extensions of, and 
comments about, it. 



A. Dielectric embedded in vacuum 

Ultimately, the physics underlying the Casimir effect 
emerges from the fact that every eigenmode of the photon 
field has zero point energy £"„ = (l/2)?iw„. The Casimir 
energy is simply the difference in zero point energies be- 
tween any two well defined physical situations 
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Sometimes we will need a regulator to make sense of this 
energy difference, though in many other cases of physi- 
cal interest (as in dielectrics) the physics of the problem 
will automatically regulate the difference for us and make 
the results finite. But adding over all eigen-modes is pro- 
hibitively difficult, so it is in general more productive to 
replace the sum over states by an integral over the density 
of states. 
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Suppose we have a finite volume V of some bulk di- 
electric in which the dispersion relation for photons is 
given by some function w(fc), which describes the pho- 
ton frequency as a function of the wave number (three- 
momentum) k. To keep infra-red divergences under con- 
trol introduce a regulator by putting the whole universe 
in a box of finite volume T4o • Then to calculate the total 
zero-point energy of the electromagnetic field we simply 
have to sum the photon energies over all momenta (and 
polarizations), using the usual and elementary density 
of states: [Volume] d^k/{2TT)^. (In the next section we 
shall look at finite-volume corrections to this density of 
states.) 



Including photon modes both inside and outside the 
dielectric body is achieved by calculating 
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Note that outside the dielectric we have simply taken the 
photon dispersion relation to be that of the Minkowski 
vacuum uj^k) = ck. 

At low frequencies, we know that the dispersion rela- 
tion for a dielectric is simply summarized by the zero- 
frequency refractive index n. That is 



Lu(k) — > ck/n 
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On the other hand, at high enough frequencies, the pho- 
tons propagate freely through the dielectric: They are 
then simply free photons traveling through the empty 
vacuum between individual atoms. Thus 
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In the absence of the dielectric, we can calculate the 
total zero-point energy of the Minkowski vacuum as 
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Let us now subtract the two zero-point energies given in 
eqns. (0) and (pi) as 
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which defines Ecasimir for a dielectric body embedded in 
vacuum. 

From the above asymptotic analysis, we know that the 
integrand must go to zero at large wave-number. In 
fact, for any real physical dielectric the integrand must 
go to zero sufficiently rapidly to make the integral con- 
verge since, after all, we are talking about a real physical 
difference in energies. To actually calculate this energy 
difference we require a suitable physical model for uj{k). 

Schwinger's calculation reported in Ref. ||l| is equiva- 
lent to choosing the particularly simple model 

uj{k) ^ — Q{K - k) + ck Q{k - K) (11) 

n 

whose physical meaning is immediately seen from the 
presence of the Heaviside step function Q{x); in addi- 
tion, K in a, wave-number (three-momentum) which char- 
acterizes the transition from dielectric-like behavior to 
vacuum-like behavior. (See figure n^.) For a dielectric 
body embedded in vacuum this particular model gives 
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FIG. 1. Schwinger's model dispersion relation. 

Note that the cut-off K describes real physics: It is a 
surrogate for all the complicated physics that would be 
required to make a detailed model for the dielectric to 
vacuum transition in the dispersion relation. 

Of course, this can all be recast in terms of a wave- 
number dependent refractive index 
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The boundary term vanishes because of the asymptotic 
behavior of w(fc), and the two substitutions k ~ uj{k)n/c 
and k — lo j c then give 
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In terms of the refractive index, of course, Schwinger's 
model for the dispersion relation becomes (see figure 0) 



i{k) = n e{K -k) + e(fc - K), 
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FIG. 2. Schwinger's model for the refractive index. Re- 
fractive index as a function of wave-number. 

and equivalently (see figure [|) 

n{u;) = n e{[Kc/n] - w) + e(cj - Kc). (17) 

We now go on to the more general case of two dielectric 
media. 
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FIG. 3. Sdiwinger's model for the refractive index. Re- 
fractive index as a function of frequency. 



B. Two dielectric media 

Suppose we now have two dielectrics to consider. We 
will deal with the situation where a dielectric body with 
dispersion relation uji{k) is embedded in an infinite slab 
of dielectric with dispersion relation uj2{k). The total 
zero-point energy is easily written down 
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If the embedded body is removed, and the hole sim- 
ply fiUed-in with the embedding medium, we can simply 
calculate the new total zero-point energy as 
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We define the Casimir energy by subtracting these two 
zero-point energies 
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The physical import of this definition is clear: The 
Casimir energy is defined as the change in the zero-point 
energy due to a change in the medium. 

From the general considerations in the previous subsec- 
tion we know that the integrand must (still) go to zero at 
large wave-number, and in fact, for any pair of real phys- 
ical dielectrics the integrand must go to zero sufficiently 
rapidly to make the integral converge. 

The integration by parts discussed above can now be 
extended to yield 
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While the difference between the refractive indices in the 
above expression goes to zero sufficiently rapidly to make 
the integral converge, it must be noted that the prefac- 
tor of oj^ implies that the net Casimir energy will be 
relatively sensitive to the high frequency behavior of the 
refractive indices. 

If we now use Schwinger's simple wave-number cut-off, 
we find that 
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Notice that depending on the precise values of the two re- 
fractive indices and the two cut-offs, we can easily change 
the sign of the Casimir energy. 

In particular, if we have a cavity containing vacuum 
embedded in an otherwise uniform dielectric ni = 1, and 
can for simplicity take n2 = n and K2 = K, in which 
case 
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Note that this is exactly the negative of the result for 
a dielectric body embedded in vacuum. [See Eq. (p2[).] 
This observation serves to drive home the fact that the 
Casimir energy makes no sense until one has carefully 
specified the two physical situations whose energy differ- 
ence is being calculated. 



C. A general cut-off 

In order to appreciate the degree to which Schwinger's 
expression is generic, notice that the dispersion relation 
can always be written in the form 



a;(/c) ~ ck + ck 
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This expression serves to define the function J-{k/K), 
and the cut-off scale K. We know from the general dis- 
cussion that T{Q) = 1, while J^(oo) = 0. It is convenient 
to normalize J- by ^(1) = 1/2, and thus implicitly fix K 
by the requirement that Ld{K) — cK{n + l)/(2n). With 
this notation K is the wave-number at which uj{k)/k has 
fallen half-way from its low-momentum dielectric value 
to its high wave-number free-space value. (See figure 0.) 




Casimir energy is proportional to this cut-off to the 
fourth power — this cut-off is real physics, and is not an 
artifact to be renormalized away. 

(5) It is only an accident of history that many, but by no 
means all, of the early Casimir effect calculations were 
for configurations where the cut-off dependence happens 
to vanish. (See ]1^.) 

(6) There is an overall dimensionless constant that arises 
from the detailed behavior of the dispersion relation as a 
function of wave-number — this remaining overall factor 
cannot be calculated without developing a specific de- 
tailed model for the dispersion. 

(7) The key physical insight underlying the Casimir effect 
is that the presence of a dielectric or conductor changes 
the photon dispersion relation and thereby changes the 
total zero-point energy. 

If we wish to tackle the problem of one dielectric em- 
bedded in another, the relevant generalizations are clear: 
we simply write down individual cut-off functions .?^(i,2) 
and cut-off scales -ft^(i,2) for each dielectric. The Casimir 
energy difference is now given by 
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Next we turn to a discussion of the finite-volume ef- 
fects. 



FIG. 4. Generic model dispersion relation. 
In terms of the wave-number dependent refractive index. 
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The Casimir energy difference is now simply written as 
1 
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This extension of Schwinger's result makes manifest 
the salient points of the Casimir energy in dielectrics: 

(1) The (dominant) contribution to the Casimir zero- 
point-energy is a bulk effect proportional to the volume 
of the displaced dielectric, 

(2) and equivalently, the bulk Casimir energy density is 
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(3) The dependence on refractive index at low momen- 
tum is explicitly known. 

(4) There is an explicit ultraviolet cut-off A', and the 



D. Finite- volume effects 

We now look at contributions arising from sub- 
dominant finite-volume corrections to the density of 
states. The key observation here is that the existence of 
finite- volume corrections proportional to the surface area 
of the dielectric is a generic result. The fact that Milton 
et al. encountered a surface-tension term proportional to 
(cut-off)"^ is hereby explained on general physical grounds 
without recourse to special function theory. 

Now, the fact that the dominant contribution to the 
Casimir energy is proportional only to volume is of 
course, merely a reflection of the fact that the canoni- 
cal bulk expression for the density of states is propor- 
tional to volume: [Volume] (i^fc/(27r)'^. It is reasonably 
well-known, though perhaps not so elementary, that the 
density of states is in general modified by finite volume 
effects. Thus in general we should write 
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These are the first two terms in an asymptotic expansion 
in l//c. A discussion of the general existence of such terms 



can be found in the textbook by Pathria [Q, while a 
more extensive discussion has been given in the paper by 
Bahan and Bloch |1^. For Dirichlet, von Neumann, and 
Robin boundary conditions the dimensionless variable ^ 
is a known function of the boundary conditions imposed. 
For dielectric junction conditions, the case of interest in 
the current problem, the situation is considerably more 
complicated and will be dealt with in a forthcoming pa- 
per. For now we just point out that, introducing separate 
quantities ^m and S,out for the density of states inside and 
outside the dielectric body, the first finite-volume contri- 
bution to the Casimir energy is of the form 
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For a simple wave-number cut-off, a la Schwingcr, this 
becomes 
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which turns out to be in full agreement with the results 
enunciated in pO|-p^. To actually complete the calcula- 
tion we need to evaluate both S.in and ^out as functions 
of refractive index for both TE and TM modes. This will 
be the subject of a forthcoming paper |g]. 



E. A potential source of confusion 

A potential source of confusion should now be noted 
before it leads to trouble: If we introduce a naive cut-off 
in frequency (energy) rather than wave-number (three- 
momentum) then the Casimir zero-point energy appears 
to have a different behavior as a function of refractive in- 
dex. This is merely a reflection of the fact that frequency 
cut-offs are in general ill-behaved, commonly leading to 
multi-branched dispersion relations. 

For definiteness, consider what at first sight would 
seem to be the frequency cut-off version of Schwinger's 
model 
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Here Q{x) is again the Heaviside step function, while 
ujq is now a frequency which characterizes the transi- 
tion from dielectric-like behavior to vacuum-like behav- 
ior. (See figure ^.) If we try to invert this function to ob- 
tain the dispersion relation uj(k), we see that it is double- 
valued in the range k € [uJo/c,nLUQ/c\. Momenta in this 
range will be doubly-counted, and lead to unphysical re- 
sults. (See figure |[) 
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FIG. 5. A naive frequency cut-off for the model disper- 
sion relation: Note the unphysical double- valued nature of 
the dispersion. 



For this particular model, we get 

1 /clIq \^ r s 1 

Ecasiniir = ^^^^^ [ ) P ^J ' 



(34) 



Note that the double-counting has been sufhcient to 
change the sign of the bulk Casimir energy! (See eq. 
(12)). Physically, the sign as determined in Schwinger's 
calculation is the correct one as can be seen from a gen- 
eral argument: Suppose only that the dispersion rela- 
tion is single- valued and that Lo{k) < ck, (equivalently 
n{k) > 1), then the sign as determined by Schwinger is 
correct. 

The [n^ — 1] behavior encountered above is in fact 
generic for frequency-based or energy-based cut-offs. In 
particular, (as we shall see below) if we impose point- 
slitting in time as our regulator, and then define en- 
ergy differences by subtracting two time-split energies 
which are time-split at the same physical time, then 
we are effectively imposing a frequency-based regulator 
and will obtain the [n^ — 1] behavior. On the other 
hand, if we time-split using "optical time" , the parame- 
ter Topticai = Tphysicai/n, then wc are effectively imposing 
a wave-number-based regulator, and we generically get a 
[(1/n) — 1] behavior. 

A wave-number based regulator, following Schwinger's 
prescription, is by far more physical than a naive 
frequency-based regulator. The only reason we belabor 
this point is because many calculations are carried out us- 
ing naive time-splitting in physical time, and should be 
slightly modified (so that we time-split in optical time) 
before being compared to Schwinger's result. Otherwise 
one is led into meaningless results. 

We emphasize that this is a minor side issue that does 
not affect issues of volume dependence versus surface de- 
pendence, and thus does not alter the cut-off dependence. 
At worst these technical issues influence the behavior as 
a function of refractive index, and even then the order 
{0[n— 1]) of the Casimir energy for dilute media (n « 1) 
will not be affected. 



III. PHYSICAL DESCRIPTION OF THE 

CALCULATIONS AND EXISTENCE OF THE 

VOLUME TERM 



The calculations presented by Milton et al. |10-lJ] 
avoid the density-of-states point of view by explicitly cal- 
culating the eigen-niodes (equivalently, Green functions) 
for a specific model configuration: a dielectric ball (of di- 
electric constant ei and permeability ^i) embedded in an 
infinite slab of (different) dielectric constant £2 and per- 
meability ^2- They then explicitly integrated over these 
eigen-modes to calculate the Casimir energy. 

Note that an important limitation of any such calcula- 
tion is that while the density-of-states argument applies 
to dielectrics of arbitrary shape, any attempt at explic- 
itly calculating eigen-modes must be restricted to sys- 
tems of extremely high symmetry — such as for example 
half-spaces, slabs, or spheres. 

The basic strategy is to start with the classical expres- 
sion for the energy (where we have assumed that e and 
/i are frequency independent constants) 
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We now promote the electric and magnetic fields to be 
operator quantities, and calculate the vacuum expecta- 
tion value 
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+ i(B(0,x).B(0,f))]. (36) 



The geometry is incorporated in the calculation both via 
the limits of integration and the boundary conditions sat- 
isfied by the fields. Since the above two-point functions 
are of course divergent, they must be rendered finite by 
some regularization prescription. Milton et al. use time- 
splitting (point-splitting in the time direction), defining 
the quantity E{t) by: 



E{r) = \ 



Geornetry 



d^x[e{E{T,x)-E{Q,x)) 

+^B{t,x)-B{0,x))]. (37) 



It should be pointed out that while time-splitting is 
a very powerful and technically useful ultra-violet reg- 
ulator it has the decided disadvantage of obscuring the 
underlying physical basis of the cut-off in dielectric me- 
dia. 

The technical aspects of the analysis carried out by 
Milton et al. reduce then to calculating these two-point 
correlation functions (Green functions) by explicitly solv- 
ing for the TE and TM modes appropriate for a spherical 
ball with dielectric boundary conditions. These Green 
functions can be written as a sum over suitable combi- 
nations of Ricatti-Bessel functions and vector spherical 
harmonics. 



To avoid unnecessary notational complications, we 
schematically rewrite the above as 



E{r) 



d^xG[^^f,]{T,x;0,x), 



(38) 



where Gi^^^] (i, x; i', x') is simply shorthand for the linear 
combination of Green functions appearing above. 

These Green functions should be calculated for three 
different geometriesjj: 

Case I: A dielectric ball of dielectric constant ei, perme- 
ability /ii, and radius R embedded in a infinite dielectric 
of different dielectric constant £2 and permeability fi2- 
(In applications to sonoluminescence, think of this as an 
air bubble of radius R in water.) 

Case II: A completely homogeneous space completely 
filled with dielectric (£2,^2)- (In applications to sonolu- 
minescence, think of this as pure water.) 
Case III: A completely homogeneous space completely 
filled with dielectric (ei,/ii). (In applications to sonolu- 
minescence, think of this as pure air.) 




FIG. 6. Case I: Dielectric ball (foreground) in dielectric 
background. 



*Notice that in Milton et al. the dielectric properties of these 
media are taken to be frequency independent, with the cut-off 
being put in via time-splitting. 




FIG. 7. Case II: Homogeneous dielectric background. 




FIG. 8. Case III: Homogeneous dielectric foreground. 

We are in complete agreement with the extant cal- 
culations and results for these three individual Green 
functions — where we do not agree, as will be shown, is 
in the way that these three Green functions are inserted 
into the computation for the Casimir energy. 

Milton et al. calculate an "energy difference", which 
we will call E surf ace, and which they define as 
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(39) 



The computation of this quantity is mathematically 
correct — unfortunately this is simply not the physically 
correct quantity to be computing, and this definition is 
equivalent to explicitly excluding by hand the dominant 
volume contribution. 

The appropriate physical quantity to compute is pi 
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(40) 



Observe that this quantity is simply the difference 
in energy between two situations: (Case I) having the 



dielectric ball present and (Case II) replacing the di- 
electric ball by the surrounding medium. This is ex- 
actly the quantity that Schwinger calculates in refer- 
ence z^/. Within the context of sonoluminescence this 
is the Casimir energy released in evolving from bubble to 
no-bubble. In contrast, the definition of Egurface above 
is not a physical difference between any two real phys- 
ical situations. Rather it is a doubly subtracted quan- 
tity which does not respect the boundary conditions for 
the problem it tries to solve. In particular, E^urface is 
not physically the same quantity as that calculated by 
Schwinger. 

If we look at the difference between the appropriate 
definition (Schwinger's) and that of Milton et al. we see 
that 



^-^ ^Casimir ^surface 
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(41) 



r<R 



This difference is now easily seen to be the missing vol- 
ume term: Remember that G/// and G// are Green func- 
tions corresponding to spaces completely filled with ho- 
mogeneous dielectrics — therefore these Green functions 
are translation invariant. (When we express these Green 
functions in terms of spherical polar coordinates it is not 
obvious that they are translation invariant, but because 
the physical dielectric in Cases II and III is translation 
invariant, the Green functions must also be translation 
invariant.) This observation permits us to pull the Green 
functions outside the integral, so that 



AE^ 



^V [Giii{t, 0; 0, 0) - Guir, 0; 0, 0)} . (42) 



This explicitly shows that the term omitted in the anal- 
yses of Milton et al. is a volume term. As we show below, 
it is in fact exactly the term required to bring that cal- 
culation into conformity with Schwinger's result. 

We shall show this by computing the difference term in 
two separate ways: (1) using translation invariance, di- 
mensional analysis, conformal symmetry, and elementary 
quantum field theory it is possible to calculate this differ- 
ence term from first principles for arbitrary geometries, 
and (2) particularizing to the case of a spherical dielectric 
ball, and working in spherical polar coordinates, we shall 
use Milton et aL's own formulae for these Green func- 
tions and give an explicit expression for this difference as 
a sum over integrals involving Ricatti-Bessel functions. 
These integrals and sums will be evaluated explicitly in 
closed form, and we shall explicitly see how the correct 
volume, cut-off, and dielectric dependence emerge from 
the calculation. 



IV. AN ELEMENTARY QUANTUM FIELD 
THEORY CALCULATION 

To evaluate the time-split energy density we apply 
some standard results from flat-space Minkowski quan- 
tum field theory. Note that in the Feynman gauge, the 
two-point function for the electromagnetic vector poten- 
tial is given by 



{T{A,{x)A^{y))) 



Vp.f 



(2^)2 {x - yf 

_Jl Tlfil' 



(2^)2c2[t(x)-t(y)]2-(.f-y)2' 

and for the Minkowski vacuum 
{E{T,x)-E{0,x))^{BiT,x)-B{0,x)) ^ ^^ 



(43) 



(2^)2 (ct)4 ■ 

(44) 



The (regulated) vacuum energy is simply [cf. eq. (37)] 

E{T)vacuum = J^;^ V J^^- (45) 

However, the effect of adding a frequency independent 
dielectric constant, can be simply mimicked by re-scaling 
the physical metric rj^^, thereby introducing what we 
shall call the "optical metric" 



7°?'™' =diag[+l/n2, -1,-1,-1]; 



y optical 



diag[-|-r 



4,-1,-1] 



(46) 



(47) 



Note that on-shell photons in the dielectric will travel 
along null curves of this optical metric. The introduc- 
tion of this optical metric is a nice technical trick for 
dealing with homogeneous and non-dispersive dielectrics 
by means of viewing them as special deformations of 
Minkowski space. 

In terms of this optical metric we can still write 



E{t) 



optical 
dielectric 
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This energy is to be interpreted as that which would 
be measured by a hypothetical observer for whom the 
optical metric would be physical. This energy must be 
translated back to physical quantities by means of the 
following equivalences: 



^physical ^ ^ optical / 1^'^ 
'^physical — "^optical ^; 
^physical — ^optical • 



(49) 
(50) 
(51) 



The effect of these translations is that the time-split en- 
ergy density (time-split in physical time) becomes 
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But, on the other hand, time-splitting in optical time 
yields 
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This now permits us to evaluate exactly the missing 
piece that Milton et al. discarded; with the convention of 
time-splitting in physical time, that term is 
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Since Milton et al. express all their results in spher- 
ical polar coordinates, as sums over an infinite number 
of Ricatti-Bessel functions, it is far from obvious that 
the rather formidable expressions encountered in those 
analyses reduce to the simple and exact result displayed 
above. In the next section of this paper we will do ex- 
actly this by invoking a long and turgid agony of special 
function theory. 

Before leaving this section, however, we wish to 
emphasize that comparing this result to the origi- 
nal Schwinger calculation requires one additional mod- 
ification: Schwinger regulated his energy calculation 
by introducing a cut-off at fixed wave-number (three- 
momentum). This is equivalent to time-splitting at fixed 
optical time. To see this, recall that frequency and wave- 
number cut-offs are in a non-dispersive dielectric related 
by 



nujo — cK, 



in which case 



'^physical ^^ ^/^O^ 

while on the other hand 
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Thus, if we regularize by time-splitting, and perform the 
subtraction at fixed optical time, we have 



A rpphysical 

fixed Taptical 



rn/ \physical 

— ^{'''optical) dielectric 1 



T-i/ \physical 

^ ['''optical ) dielectric 2 



(27r) 



2 * physical 



1 1 



12 



\C '^optical ) 



(58) 



10 



This is exactly Schwinger's result as enunciated in [Q . 
(More precisely, it is as close as one is ever going to get 
considering that a sharp cut-ofF is just not the same as 
time-splitting: There is no reason to expect the dimen- 
sionless coefficients to be the same for the two calcula- 
tions. We will comment, near the end of the next section, 
how we can switch from a time splitting cut-off to a sharp 
cut-off in wave number and obtain exactly Schwinger's 
coefficient.) 

In summary, what we have shown achieved so far is to 
show the following: 

(1) Schwinger's calculation ||l| showing that there is a 
bulk Casimir energy in a generic dielectric is correct. 

(2) The analyses of Milton et al. pO|-[l4| omit the volume 
contribution to the Casimir energy. 

(3) The volume contribution that was omitted in those 
analyses can be calculated exactly, without resorting to 
special function theory or asymptotic analysis, and the 
omitted term exactly reproduces the volume contribution 
as evaluated by Schwinger. 

(4) The volume term in the calculation was missed be- 
cause of misidentification of a Green function, due to the 
incorrect application of boundary conditions. 

(5) This argument is not limited to systems of spherical 
symmetry. From the preceding it is clear that the argu- 
ment continues to hold for arbitrarily shaped dielectrics. 

In the next section we will evaluate /S.E given in eq. 
( [4l| ) in yet a different way: by working directly from 
Milton's expression. 



V. THE ENERGY DIFFERENCE TERM IN 
SPHERICAL COORDINATES 

From the preceding discussion we have isolated the 
term in the Casimir energy that Milton et al. 10-1^] 
omitted. It is precisely the difference between the Green 
functions (of Case III and Case II) integrated over the 
interior of the bubble (See eq. (p])). We can thus imme- 
diately jump into the middle of the technical computation 
and directly evaluate this difference term. This calcula- 
tion has the advantage of quickly getting to the heart of 
the matter. 

From the two recent papers ||l^ , eq. (41)|j, or from [|2| , 
eq. (4.2b), we can write 
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X Xiii{kui-,r) - Xn{kii,r) 



(59) 



^Note that this is eq. (42) in the hep-th version. 



Here we use the notation k = |w|n, with n the appropri- 
ate refractive index, and define the quantity Xjj{kii,r) 
by 



XH{kH,r) = ^(2i'+ 1)<^ 2 k]j Fl\kii-r, 
£=1 [ 

J.2 Qj. 



-4j^'[FPikir.ry)] 



dr 

[{Fn - (Gf )] 



(60) 



The Ff and Gi Green functions are identical to those 
in |[10|-p2[. Notice that when making the substitution 
(//) — > {III) we should also change the refractive index 
that implicitly appears in the factor k. 

In writing these equations we have used the fact that 
the energy difference is known to be a real quantity, so 
we can immediately discard the imaginary part of the 
above expression without bothering to explicitly evaluate 
it. There is no subtle physics involved here — this is just 
the standard procedure of using complex exponentials to 
represent a harmonic time dependence and then taking 
the real part at the end of the calculation. In particular 
there is no physics hiding in the imaginary part of the 
above expression. Not only is there no need to calculate 
the imaginary part, but it is physically meaningless to do 
so. 

If we were integrating over all space, then — follow- 
ing |10|-[l^ — the partial derivative terms in X{k,r) 
could be safely dropped. Because we are only integrating 
over a finite region we must explicitly keep these deriva- 
tive terms. 

Now because k — \uj\n, while Fg and Gg are functions 
of fc, they are therefore functions only of the absolute 
value of uj. This permits us to write 
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(61) 



and finally 



AE = I — cos{ujt) 



2n 



^dr 



xlm 



Xiii{kui,r) - Xn{kii,r) 



(62) 



The Green functions for Cases II and III are very sim- 
ple (since we are dealing with homogeneous spaces) and 
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the imaginary parts can easily be read off from |L0| eq. 
(29), eq. (35), or (l| eq. (3.7)| 



Im{^a = Im{Ge} = Im{tk n{k,r) h,{kr)} 
= k jjikr) 

[krY ■ 



(63) 



Here we have introduced the Ricatti-Bessel function 
si{x). (We refer the reader to the Appendix for details.) 
Inspection of the derivative pieces yields 



Im 



d_ 
dr 



d 



r—r'[F,{k-r,r')] 



dr 



= Im<! — 
or 



r—r'[G,{k;ry)] 
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= Im{ifc [{xhi{x))'\kr {xji{x))'\kr 

+ {xht{x))\kr {xjl{x))"\kr]] 

= k Re{e^(a:)|fe,. s'^{x)\kr + et{x)\kr s'f{x)\kr} 

= k {s'^{x)\kr s'^{x)\kr + Sl{x)\kr s'p{x)\kr} , (64) 

with ef,{x) another Ricatti-Bessel function. 

To derive these results we can either use brute force 
starting from eq. (35) of pl|, or alternatively we could 
inspect equations (19-22) of |10|. (Note that there was an 
overall change in normalization between the 1980 paper 
and the 1995 and 1996 papers, and take account of the 
fact that we are here retaining the total derivative term.) 

All together, this implies 

Im Xiii{kin,r)~ Xii{kii,r) 

X {2{stix)f + i[se{x)yf + Si{x)s';{x)} |fc,„, 
-[(///)-(//)]. (65) 

The sum over Ricatti-Bessel functions is now easily 
and exactly evaluated using the results in the Appendix. 
We simply get 



lm{Xiii{kiii, r) - Xii{kii, r)} = 2 



- 2{kr) 



{III) 
{II) 



= A[khi-kh] 



4 [n? - 71^] |c^|3. (66) 



Note: The two lines given in those equations are individ- 
ually the Green functions for Cases II and III (defined over 
the whole space). The combination given in ||l[J-[L2| is not a 
Green function of any differential operator. See the discussion 
surrounding eq. (^l|) earlier in this paper. 



Notice that the above expression is independent of r — a 
result that is by no means obvious from the original def- 
inition. (But we knew from the fact that the underlying 
Green functions are translation invariant that this quan- 
tity had to be independent of position at the end of the 
day!) 

Turning again to the energy difference 



AE = [nl -nl] — [cos(wt)] / r^dr 

Jo 27r Jq 
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Which is now easily evaluated, using x — cut, as 
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X cos{x)dx. 



As a penultimate step we use 

r'+oo 



X cos{x)dx = 6, 







and V = (47r/3)i?^, to write 



AE = ^[nl-r.l]l,V. 



(68) 



(69) 



(70) 



Re-inserting factors of h and c, which have been sup- 
pressed for clarity, we see 



AE{Tphysicai) ^ — he [nl - nl] 
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physical 
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Note that this is exactly the same result {ineluding nu- 
merieal eoeffieients) as the previously calculated from 
general considerations using translation invariance and 
the strategy of introducing the optical metric. [See eq. 

Finally, as discussed previously, we note that to com- 
pare this result to Schwinger's we should be time-splitting 
in optical time, and so must absorb a few factors of re- 
fractive index into the time-splitting parameter. [See eq. 
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This is exactly the time-split version of Schwinger's ex- 
pression (^4|) for the bulk contribution to the Casimir 
energy. 

We have thus evaluated the difference term, AE = 
Ecasimir — E surf ace, in two independent and completely 
independent ways. The two calculations agree down to 
exact numerical coefficients. Turning this result around, 
we have 



E, 



Casimir 



= AE + E. 



surface- 



(73) 



So we see that the term we have just calculated provides 
the bulk volume contribution to the Casimir energy that 
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is required to bring Milton et al.'s calculations into con- 
formity with Schwinger's arguments ||l|, and into confor- 
mity with the general arguments provided in this paper. 
To finally compare the absolute normalization of this 
result with Schwinger's, we must explicitly replace the 
time-splitting by a wave-number cut-ofF. Simply go back 
to eq. d63) , introduce a wave- number cut-off, and let the 
time-splitting parameter go to zero, to obtain 



A^ 
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— / r''dr 
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xlm{f{kiii)Xiii{kiii,r) - f{kii)Xii{kii,r)}. 

(74) 

The function f{k) is any real smooth function with 
/(O) = 1 and /(-|-oo) = 0. For simplicity we have tem- 
porarily used the same cut-off function for the two media. 
We generalize the argument below. The evaluation of the 
various pieces of the Green functions proceeds as before, 
so that eq. now mTi) becomes 



AE= — 



2^ 



r^drA{f{kiu)k]ij - f{kii)k^j}. 



(75) 



After a change of variables, and an explicit evaluation of 
the volume integral, we easily get 



AE 



1 



1 

n2 



V 



d-^k 
(2^ 



f{k)k. 



(76) 



Approximating air by vacuum, (i.e. setting ni — 1), and 
inserting a sharp cut-off at wave-number K, (by setting 
f{k) — Q(K—k)), this is exactly Schwinger's bulk volume 
term for a cavity in a dielectric [eq. (pq)] down to the 
last numerical prefactor. 

Now the use of a single cut-off function, while mathe- 
matically more transparent, is physically dubious. If we 
introduce separate cut-off functions for the two media the 
modifications are straightforward. First 
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(77) 

The functions fa{k) are now any two real smooth func- 
tions with /a(0) = 1 and /(i(+oo) = 0. We now have 

/>+oc 7 pR 

AE^ -^ / r^drA{fin{kni)khi-fH{kii)k],}. 
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(78) 

After a change of variables, setting //// — /i and /// = 
/2, and explicit evaluation of the volume integral, we 
easily get 



AE = V 



d^k 

'{2^ 



h{k) hik) 



ni 



"-2 



k. 



(79) 



Inserting sharp cut-offs in wave-number, a la Schwinger, 
is accomplished by setting 

fa{k) = [OiKa - fc) + n,e(fc - Ka)] e(ifoo - k). (80) 

These cut-off functions may look a little mysterious, de- 
pending as they do on three cut-off scales Ki, K2, and 
Koo- The two cut-offs A'l and K2 are physical — they 
describe the wave-numbers at which the two dispersion 
relations effectively approach the vacuum dispersion re- 
lation. The third cut-off K^o is purely a mathematical 
convenience to make the integrals converge. K^o should 
be taken to be much larger than cither Ki or K2- With 
this notation 
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(81) 



That is to say, these particular cut-off functions have 
been chosen to mimic the model dispersion relation used 
in the mode sum discussion (section H) . 

Evaluating the integrals, and reinserting factors of Ti 
and c, we can write 
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(82) 



Note in particular that K^o has quietly disappeared from 
this answer. This is exactly the same as the result ob- 
tained by mode sum arguments in section ^ down to the 
last numerical prefactor. [See eq. (p^).] 

A slightly more general cut-off is to use the wave- 
number dependent refractive index to construct the cut- 
off according to the scheme 
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(83) 



In which case, recalling that that Ua = na{k = 0) 
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(84) 



The Casimir energy difference, including reinserted fac- 
tors of fi and c, is now 



AE^hV 



d^k 

1^ 



[ui{k)-uj2{k)]e{Koo-k). (85) 



We can now let Krx, tend to infinity, and thereby recover 
eq. (|0|). 

In summary, the energy difference term calculated via 
explicit summation over Ricatti-Bessel functions is com- 
pletely in agreement with the dominant volume term 
coming from simple mode sum and density of states ar- 
guments. 
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VI. FULL CALCULATION FROM FIRST 
PRINCIPLES 

A. Background 

There is also a certain amount of value to going to the 
trouble of redoing the calculation from first principles. 

Recall that the physical quantity that we are interested 
in evaluating for the calculation of the Casimir energy is 
the total energy. The total classical energy is computed 
from the integral of the energy density over the geome- 
tryQ 



E 



<j Geornetrv Z \ 



(86) 



After quantization this energy is related to the vacuum 
expectation value of the squares of the electric and mag- 
netic fields in the dielectric. 

Since the electric and magnetic fields are related to 
derivatives of the electromagnetic potential, the energy 
density is calculable in terms of the two-point correlation 
function of the electromagnetic field; these correlators are 
in general divergent, and we need to perform renormaliza- 
tions before comparing with the results of experiments. 
This, of course, is nothing but the general recipe for eval- 
uating quantities in any quantum field theory. 

This calculation has been done in detail by Milton et 
al. pO|-p^ for the geometrical configuration previously 
described. They computed the total energy E by inte- 
grating the electromagnetic energy density over the ge- 
ometry and then rendered the expressions finite using 
point-splitting regularization in the time direction. We 
quote the expressions from them (see below for the vari- 
ous definitions), as 



S = Rc 
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(87) 



Here the total derivative can be shown (Ref. 10 of Mil- 
ton) to integrate to zeroJ^ In this equation the w-integral 



'^This expression is valid for frequency independent e and /x 
only. 

^We reiterate that this vanishing of the total derivative terms 
depends crucially on the continuity and differentiability of the 
Green functions. 



is for the regulator in the time-direction, and the ra- 
dial functions Fe{r,r') and Gi{r,r') are related to the 
Green functions for the electric and magnetic fields in 
the appropriate geometry. For a dielectric sphere of ra- 
dius R, permeability ei and permittivity /ii, embedded 
in another dielectric with 62 and /i2, using the results of 
Boyer [|l^i|8| as adapted by Milton et al., it is possible 
to readily check that these objects are the radial part of 
the Green functions for the wave equation in spherical 
coordinates. 



+ k^)F,{k;r,r') 
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and are given explicitly as follows: 
For r, r' < R, 

Fi{kin; r,r')= ikinjt{kinr<) 

x[/i^(fci„r>) - Apjtihnry)]. 

For r, r' > R, 

Fi{kout;r,r') = ikouthi{koutr>) 

'x[ji{koutr<) - Bphi{koutr<)]. 



(89) 



(90) 



A similar expression holds for Gi. In these two ex- 
pressions the wave numbers are kin = |w|^ei/ii and 
kout = \^\\/^2fJ-2- The function ji{x) is the spherical 

Bessel function of order i and hi{x) = h^ '(x) is the 
spherical Hankel function of the first kind. The quanti- 
ties Ap, Aq, Bp and Bq, axe given in |10|-p^ and their 
explicit forms are obtained by requiring that the Green 
functions be solutions to the appropriate equations sat- 
isfying the correct boundary conditions. 

Milton et al. in their papers give the expression for the 
energy that we quote above, eq. (|87|), after they have 
integrated over all of space before subtracting the ener- 
gies for each material configuration. This is dangerous 
since there is then an additional infrared divergence aris- 
ing from the volume integral which is not regulated by 
the time splitting-procedure, which only regulates UV 
divergences. To avoid these unnecessary pitfalls we will 
write down the energy density, perform the subtractions 
at the level of the densities (local subtractions) and then 
integrate. 

Using the spherical Bessel function form for the electric 
and magnetic fields satisfying the Maxwell equations, the 
energy density, T", is seen to be given by 



r"(r) = Re 



-|-oc 



2^' 



Xik,r] 



(91) 



where we have used the notation k — \uj\n with n the 
refractive index appropriate for the medium, and defined 
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the quantity X{k,r) via 



X(fc,r) = ^(2£+1)< 



e=i 



,2 , ^(^+1) 



[Fi{k;r,r)] 






ri=r2=r 



(92) 



This is of course equivalent to eq. (|87|) above. Here 
ki = \uj\ni represents the wave-number for Ught of fre- 
quency a; in a medium with index of refraction n/. Sim- 
ilar results, with appropriate substitutions for the mo- 
menta and Green functions, hold for the situations when 
we make reference to Cases II and III. Note that when 
making the substitutions (/) -^ {II) or (/) -^ (HI) we 
must also change the refractive index that implicitly ap- 
pears in the factor k. In addition, it should be borne in 
mind that fc/ is a function of position: kj — ni|a;| = kju 
inside the dielectric sphere, whereas kj — n2\w\ — kjj 
outside the dielectric sphere. 

The Casimir energy is now obtained by taking the dif- 
ference in energy densities and integrating over all of 
space while paying attention to the appropriate index 
of refraction for each region of space. For example, the 
Casimir energy difference between Cases I and II is given 
by 



^Casimir -tle'^ 



r dr —— 



2^ 



xe-''^^X{ki,kn,r) 

where we have defined the quantity X{ki, fc//, r) by 
X{ki,kii,r) = Xi{ki,r) - Xii{kii,r) 



(93) 



^J2i2i+1)< 



kf + 



X [F^{ki;r,r)+Gl{ki;r,r)] 

Id d r^r,, X 

H — n -7^—riT—r2\Fi (fc/;ri,r2) 
r"^ OTi or2 



+Gl{kj;rur2)] 
m - {II)] ■ 



r-i=r-2=r 



(94) 



To turn this into eq. (41) of reference ||ll[| requires a few 
technicalities which were not made explicit by Milton but 



we will show here. (Note that the above is his eq. (42) 
in the hep-th version of the paper.) 
Observe that 



d_ 
dr 



d 



r^r' [Fl{ki-r,r') +G\{ki-r,r')\ 



dr 

d d 
^^—ri^T—r2 [Fg{ki;ri,r2) + G^(fc/; ri, r2)] 
ori or2 

+r-^r'[Fl(ki-r,r') + Gl{ki-r,r')\ 

d d 
-7rrri—r2 [F/(fc/; ri, ra) -I- G^(fc/; ri,r2)] 



rx=ri=r 



t{i+l) 



r^[Fl{ki-r,r)+Gi{ki;r,r)\ 



-25(0), 



(95) 



where in the last line we have used the differential eq. 
( psi ) satisfied by the Green functions to replace the double 
derivative at the same point by an explicit polynomial. 
The term proportional to 5{Q) arises from the source 

defining 



term on the RHS of the differential equation 
the Green functions. Fortunately this term is completely 
independent of the dielectric properties and in fact is in- 
dependent of all properties of the medium. Thus when 
calculating energy differences this term cancels identi- 
cally. (Earlier calculations often quietly discard this term 
without even mentioning its existence.) 

Inserting this into the general formula for the Casimir 
energy, eq. (|9^), yields eq. (41) of reference |11|. Explic- 
itly, it allows us to re-write X{ki, kji,r) as 

X(fc/,fc7/,r) = ^(2£-M)i2fc2[F/(fc7;r,r) 



+Gliki;r,r) 



1_ d_ 

J.2 Qj, 






F^iki;r,r') 



+Gi(fc/;r,r') 
lil) - (11)] , 



(96) 



which is now straightforward to evaluate. 

We know that k is only dependent on the absolute value 
of Lu, which permits us to write 



r'+oc 



duj 



Ecas^m^r = Re<( — / r^ dr / ^ [g-^^^ ^ -+'"^1 



®The quantity X{k,r) defined here is (apart from a phys- 
ically irrelevant 5(0) term, see below) equal to the quan- 
tity X{k,r) defined in the previous section — these quantities 
differ only by the way that the derivative terms have been 
manipulated. 



xX{ki,kH,r) 



Wherefore 

^Casimir — ^^ 



(97) 



/ ■ r 2, r°° "^ < ^ 

< —I \ r dr — cosIlut) 

[Jo Jo 2^ 
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xX{ki,kii,r) 



(98) 



We are thus led to 



E, 



Casimir 



+ OC 



Jo 27r 
X cos(cjt) lm{X{ki, kji, r)}. (99) 



With due caution^ the relevant pieces of the Green func- 
tions can be read off from [IT0HT2I . 



B. Calculation of the various Green functions. 

The general form of the Green functions were given 
above, but one must require that the boundary condi- 
tions appropriate to the geometry and the physics be 
correctly incorporated into them. This means that they 
must satisfy appropriate continuity conditions derived 
from Maxwell's equations. In terms of the fields, 



E±, eEr, —B±, and B^, 



(100) 



must be continuous. In terms of the radial Green func- 
tions, Fi and Gi, 



d 

^iFi, Gi, -^rFi 
or 



and -—-rGi 
e or 



(101) 



must be continuous. 

We now deal with each of the above three cases. 
Case I: 

Consider a spherically symmetric bubble of radius R and 
index of refraction ni embedded in an otherwise homo- 
geneous material of index of refraction n2. (See figure 
y). The Green functions satisfying the correct boundary 
conditions are: 
For ri, r2 < R, 



F/, Gi{ri,r2) = ikiii jt{kiiir<^] 
X [htikiiir^) - 

For ri,r2 > R, 



A'f^g 3i{kiiir>)\ , (102) 



F' ,G[[ri,r2) ^ikii hi{kiiry) 

X [jiikiir^] - B%a h,{knr<)\ . (103) 

(See equations (12a) and (12b) of |jlffl, and eq. (16) 

of fl^ .) The quantities Ap q and Bp q are the ones given 

by Milton et al. |l^-|l|. 

Case II: 

For the configuration that we have called Case II (see 

figure 0), the Green function is: 

For all ri , r2 

F/^Gf (ri,r2) = ikn jdkiir<) he{knr>). (104) 



Case III: 

Finally for Case III (see figure ||) the Green function is: 
For all ri , r2 

FP', Gf ^(ri, r2) = ikju Je{kiiir<) h,{kiiir>). (105) 

We are now ready to explicitly compute the Casimir en- 
ergy 

The imaginary parts can now easily be read off. For 
Case I with r < R, some straightforward algebraic ma- 
nipulations lead to: 

Im{F/}(fc„/;r,r) = Im{G,^}(A:///; r, r) 



= Rclkjjj je{kjjjr) 

X [hi(kiiir) - Ap^Q ji{kiiir)] | 

= (1 - RciA'pa}) km j^{kjur) 
sjikiiir) 



(1 - Re{A'pa}) km 



{kmrf 



(106) 



where to simplify the writing we have again introduced 
the Ricatti-Bessel functions. (See the Appendix.) 
Similarly, for r > R we get 

Im{F/}(fc//; r, r) = Im{G^^}(A://; r, r) 
= Relkii hi{kiir) 

X [jiikur) - B'pq he{kiir)] | 
= kii jUkiir) - kii Re {s|, ^ hj{kiir)} 



^11 



{kiirY 



kii 



^^{B'F.G^ikiir)] 



{kiir) 



(107) 



(108) 



On the other hand, for Case II we simply have 

Im{F/^}(A://;r,r) = Im{Gf (fc//; r, r)} 

= Re{kii jt{kiir) hi[kiir)} 
= kii ilikiir) 

sjikiir) 

" kii 71 T2"- 

(kiiry 

Inspection of the derivative pieces yields, for r < R 



Dl^lm\—r^—r2[F,pr^,r2) 

[or I Or2 - ri=r2=r 

= lm{ikm {xji{x)y\kj,jr 

X [{xh({x)y\krjjr-A'pG ixje{x)y\krjir]} 

= km Re{s'g{x)\k,iir 

X [Sl{x)\kuir — ApQ Sg{x)\kijir\} 

= (1 - RelA'p^c}) km {Weix)]\,,,r} (109) 
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While for r > R we have: 

Dint = ^"^|a;:^^ia^^2 [i^/(ri,r2)] 

= lin{ikii {xhi{x))'\kjir 

X [{Xjiix)y\kj,r-B'pa {xhi{x)y\kj,r]} 

= kii Re{ei{x)\kur 

^kii {4(a;)|fc„r Si,{x)\kj,r} 

-kii Yle{Bp(. e'i{x)\kiire'i{x)\kiir}- (HO) 

Finally, for Case II we have the relatively simple result, 
valid for all r, that: 

= Im{z/c// {xhi{x))'\ki,r {xji{x))'\ki,r} 

= A:// Re{e^(a;)|fe„r 4(2;)U„r} 

= A:,, {[4(x)]2|fe„,,}. (Ill) 



n^^ = Tm 

all space — 



When calculating differences, a large number of pieces 
in these expressions quietly cancel. What remains is still 
complicated, but using some identities between Bessel 
functions and their derivatives such as 



2 *'^(^) = ■5"(^) + ^t{x)-> 



(112) 



and an identical equation that holds for ee,{x), the calcu- 
lation simplifies considerably and gives 

Qin = ln\{Xi{ki,r)}in 

e.=\ 
X {l\si{x)f + \s\{x)f + si(x:)s",{x)] |fc,„,. (113) 

Rearranging this yields 

Qin = lui{Xi{ki,r)}in 

f=i 
X {2[se{x)]^ + [4(x)]2 + si{x)s';{x)} !,.,„, 

, CXD 

-^j:(^e + me{A'p + A'a} 

i=\ 

X {2\st{x)f + \s',{x)f + s,(a:)4'(x)} |fc,„,. 

(114) 

The first two lines are immediately recognizable as the 
terms we encountered previously while calculating the 



energy difference, Ai? = Ecasimir - E surf ace- For these 
two terms the sum over £ can be performed exactly with 
the result that 

Qtn = Im{X(fc/, /c//, r)}in 



= 4fc3 



km 



Y^{2e+l)Re{A<p + A'a} 
1=1 

X {2[,S,(.t)]2 + [S',{X)]^ + S,{X)S';{X)} \k,,,r. 

(115) 

Note that for the Ap q terms we cannot explicitly per- 
form the £ summation because of an implicit (, depen- 
dence in the coefficients A^pQ. Also notice at this point 
that the A\^ g-terms in the above expression were the 
only pieces retained in ||ll| , [l2| . 

Taking the cue from the above result, for the correct 
result in the region inside the bubble then we can write 
(using self explanatory notation) 



= 4n3 |c^|3 + Q|„"'■>=^ 



(116) 



Turning next to the region outside the dielectric 
sphere, we have 

Qout = Im{X/ (fc/,r)}o„t 

1=1 

X {2[si{x)Y + [3'i{x)f + St{x)4{x)} \k,,r 
J oc 

-^Y.(^e+me[{B'p + B'^) 

£=1 

X {2[e,{x)f + [e',{x)]^ + eiix)e'l{x)} ] U„. 

, OO 

kjj 
r2 



-^kh-^Y.(^i+mem + B'a) 

e=i 

X {2[ee{x)f + [e',{x)f + e,{x)e'l{x)} ] U„.. (117) 



That is 






(118) 



Now looking at the Green function for Case II, we get 
the simple result that 

, OO 

Im{Xu{ku,r)}ali space = 2-^ ^(2£ + 1) 



X {2[si{x)f + [s'e{x)f + si{x)4{x)} |fe„^ 
// — -±"-21^1 • (119) 



= 4fc?, = 4n?|cjP 



Turning again to the total Casimir energy, we can eas- 
ily perform the subtractions and substitutions, to obtain 
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E, 



Casimir 














2tt 



E, 



Casimir 



+ 00 



duj 



r'dr I ""^^r./,.._.^,.3 ,„.„^,.3 



2tt 



Hf{kiii)khi - I{kn)k],} 



- / r^dr / — [cos(wt)] Q 

Jo Ja 27r 

■ / "^^dr / -— [cos(wt)] Q 

7fl Jo 2^ 



surface 



surface 
out 



(120) 

The remaining r and w integrals for the uj'^ term are easily 
performed. Reinserting appropriate factors of h and c we 
are led to: 



i-R n + oo 1 

Jo Jo ^TT 



surface 



surface 
out 



E, 



Casimir 



I H - "2] ^^(^^ 



r^dr 



+00 



duj 

2^ 



cos(wt)] hQ'l'^''^'""'{kiii,r) 



+00 



duj 



+ 1 r'^dr I '^[coB{ujT)]hQout{kii,r). 
Ir Jo ^tt 



(124) 

With suitable changes of variable, performing the r 
integration for the k^ terms, and reinserting h and c as 
appropriate, we have 

Ecasiniir = +2^ / TTTTTT o ^/(^){^1 (^) " ^2(^)1 



nf{n2\uj\)Q 



(27r)3 2 

Jo 27r 

00 r+00 I 

out ' 



(121) 



r dr 



fl 







27r 



(125) 



This is now explicitly of the form: 

(Schwinger's bulk term) + (surface contribution) + • • • 

In a certain sense this terminates our calculation, since 
Milton et al. have already calculated the object we called 
E surf ace, and explicitly shown it to be a surface term 
(plus even higher-order corrections). For dilute media, 
m « 1 « 712, Milton et al. obtain in eq. (51) of ||ri| and 
eq. (7.5) of 11 



E surf ace = "("-l " n-2)^?ic 



i?2 



+ 



1 



{cTphysicalY ' 64i? 



+ 



which is the central result of this paper. 

If we finally approximate air by vacuum, (set ni = 1), 
and make the above a hard cut-off at wave-number K 
then the first term is exactly Schwinger's bulk volume 
term, while the remaining terms are subdominant surface 
and higher-order contributions. 

With little additional trouble one can introduce sep- 
arate wave-number cut-offs for the two media in which 
case we can write 

Ecas^m^r = +2F / ^—--n{f^(k)uj^{k) - /a (fc)cc>2 (fc)} 



(122) 



(27r)3 2 

+ / r^dr / nh{nM)Q 

Jo Jo 27r 



surface 



The existence and qualitative features of this surface 
term are in complete agreement with the general anal- 
ysis adduced in this paper. 

(See also eq. (51) of Q where the related result for 
the pressure difference is presented.) 

To see what happens for a wave-number cut-off, back- 
track to eq. (99), insert a wave-number cut-off /(fc) and 
let the time-splitting parameter go to zero. Then 



r dr 



+ OC 



dui 
2^ 



n f2{n2\uj\)QlZ^ "''''' , 



(126) 



E, 



Casimir 



Jo 27r 

X lTn{f{ki)Xi{ki) - f{kii)Xii{kii,r)}. 

(123) 



Here again f{k) is any smooth real function with 
/(O) = 1 and /(-|-oo) = 0. All the computations of 
the Gree n functions remain unaltered, and we replace 
eq. (^ by 



VII. DISCUSSION 

After this relatively turgid mass of technical manipu- 
lations, the main results of this paper can be succinctly 
stated: 

(A) In a dielectric medium of dielectric constant n the 
Casimir energy is, in the bulk medium, dominated by a 
volume term: 



e: 



'hulk 
Casimir 






1 
1 

n 



(127) 



This result is completely in agreement with Schwinger's 
argument in |l|, and can be adduced from simple mode 
sum arguments as presented in section ||. We have 
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checked these mode sum arguments against general field 
theoretic arguments in sections OI and |^. 
(B) In addition, we pointed out in section y that there 
will be a sub-dominant contribution to the Casimir en- 
ergy that is proportional to the surface area of the dielec- 
tric. This surface contribution takes the generic form 



e: 



surface 

Casimir 



67r2 



ShcK'^ 



fo 



'^out 



(128) 



This term will be sub-dominant provided the size of the 
dielectric is large compared to the cut-off wavelength. 
More specifically, provided 



V/S > 1/K = \o/{2tt). 



(129) 



(C) In general, following section ||, we can expect these 
to be the first two terms of a more general expansion 
that includes terms proportional to various geometrical 
invariants of the body. By analogy with the situation 
for non-dispersive Dirichlet, von Neumann, and Robin 
boundary conditions p5| ], we expect the next term to 
be proportional to the trace of the extrinsic curvature 
integrated over the surface of the body. 

(D) These very general considerations are buttressed by 
an explicit re-assessment, in sections ^ and VI, of cur- 
rently extant calculations for a dielectric sphere. We find 
that some calculations have used an inappropriate sub- 
traction scheme to define what is taken to be the Casimir 
energy. When this is fixed, all calculations fall completely 
into line with the general considerations presented in this 
paper. Furthermore, the corrected calculations are then 
also seen to be in complete agreement with Schwinger [nl . 

(E) In previous analyses of the Casimir effect it has been 
common to fixate on the van der Waal's forces as the un- 
derlying physics ultimately responsible for the Casimir 
effect. We strongly disagree with this point of view, 
and emphasize, rather, that the presence of the dielec- 
tric medium induces a change in the dispersion relation 
and a change in the density of states, that results in a 
change in the total zero point energy. This change in the 
total zero point energy is the Casimir energy. 

For a finite-volume dielectric, relative to vacuum. 



Ea 



astmir 



2V 



(fk 1 

^2^2 



Ti [Lo{k) — ck] 



(130) 



where the dots represent terms arising from higher-order 
distortions of the density of states due to finite- volume 
effects. 

(F) We close with what is perhaps a trivial point that we 
nevertheless feel should be made explicit: The volume 
contribution to the Casimir energy is always there, and 
always physical, but sometimes because of the specific 
physics of the problem, it is safe to neglect it. 

On the one hand: Suppose we are provided with a 
fixed number of dielectric bodies of fixed shape (in par- 
ticular, of fixed volume), and suppose that we simply 
wish to move the bodies around in space with respect 



to each other. Then the bulk volume contributions to 
the Casimir energy, while still definitely present, are con- 
stants independent of the relative physical location of the 
dielectric bodies, and so merely provide a constant off- 
set to the total Casimir energy. If all we are interested 
in is the energy differences between different spatial con- 
figurations of the same bodies then the various volume 
contributions can be quietly neglected. 

On the other hand, an equally physical situation is 
this: The volume contribution is of critical importance 
whenever we want to calculate the energy difference be- 
tween an inhomogeneous dielectric and a homogeneous 
dielectric wherein the irregularities have been filled in. 
This is the physical situation for example in the case of 
bubble formation in a dielectric medium. 
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APPENDIX A: RICATTI-BESSEL FUNCTIONS 

This appendix collects some useful identities involving 
sums of Ricatti-Bessel functions. We start by noting the 
very useful result given in Abramowitz and Stegun, page 
440, eq. 10.1.45. In terms of spherical Bessel functions 



J2i2£ + i)jii^)My) = 



e=o 



sin(a; — y) 

(x-y) 



(Al) 



An equivalent result, in terms of the ordinary Bessel 
functions, can be obtained by taking the real part of 
equations 8.533.1 and 8.533.2 on page 980 of Gradshteyn 
and Ryzhik. 

For the purposes of this paper, it is more convenient 
to use the Ricatti-Bessel functions defined by 



Si{x) = xje{x) 



ei{x) — xh{(x) 



■h+i/2{x)- 



H. 



f+l/2 



(x). 



(A2) 



(A3) 



Furthermore, since the sums occurring in the current 
problem always run from £ = I, rather than ^ = 0, it is 
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useful to pull the H - 
thus write 

oo 

^(2£+l)s,(x)s,(y) 



terms to the right hand side and 



xy- 



y) 



(x-y) 



sin(a;) sin(y). 



These last two equations are the only sums of Ricatti- 
Bessel functions we will actually need. By the manner 
in which we have obtained them we can easily see that 
they are simple generalizations of the usual summation 
theorems for ordinary Bessel functions. 



(A4) 



By taking derivatives with respect to x and y it is now 
straightforward to show that 

X](2^ + l)s',ix)s',{y) = xy '^i-^f + ''°'('' " y^ 

fcosfa; — u) sinfx — -y) 
+2xy V «^ \ yj 



[x - yY {x - yY 

cos(a;) cos(j/), (A5) 



and 



E/„ „ , N / N // / X sin(x — y) 

{2t+l)s,{x)4{y) = -xy ^ ^^ 



-2x^ 



{x-y) 

cos(a:: — y) sva{x — y) 



{x-yY {x-yY 

+ sin(x) sin(y) . ( A6) 

Taking coincidence limits (y -^ x) we now get 

oo 

^(2£ + \)sii{x)sii{x) =x^ - suY{x), (A7) 



[1] J. Schwinger, Proc. Nat. Acad. Sci. 90, 2105-2106 (1993). 

[2] J. Schwinger, Proc. Nat. Acad. Sci. 90, 4505-4507 (1993). 

[3] J. Schwinger, Proc. Nat. Acad. Sci. 90, 7285-7287 (1993). 

[4] J. Schwinger, Proc. Nat. Acad. Sci. 91, 6473-6475 (1994); 

[5] J. Schwinger, Proc. Nat. Acad. Sci. 89, 4091-4903 (1992); 

[6] J. Schwinger, Proc. Nat. Acad. Sci. 89, 11118-11120 
(1992); 

[7] J. Schwinger, Proc. Nat. Acad. Sci. 90, 958-959 (1993); 

[8] In preparation. 

[9] S. Weinberg, Rev. Mod. Phys. 61, 1-23 (1989); 
[10] K. Milton, Annals of Physics (NY), 127, 46-91 (1980). 
[11] K. Milton, Casimir energy for a spherical cavity in a 
dielectric: toward a model for Sonoluminescence?, in 
Quantum field theory under the influence of external 
conditions, edited by M. Bordag, (Tuebner Verla gsge- 
sellschaft, S tuttgart, 1996), pages 13-23. See also hep 



th/951009l[ Warning: the equation numbers (though not 



£=1 



while 



oo _. 

J2i2£+lH{x)s',{x)^-x'+siiY{x 



[12] 
[13] 



the text) differ between the published version and the 
hep-th version — all equation numbers after (32) are in- 
cremented by 1 in the hep-th version. 
K. A. Milton and Y. J. Ng, Casimir energy for a spherical 
cavity in a diele ctric: Applications to Sonoluminescence, 



hep-th/9607186 



(A8) 



i=i 



and 



oo _. 

^(2£ + l)si{x)s'[{x) = - o^' + sin(a;)2. (A9) 



and 



C. E. Carlson, C. Molina-Pan's, J. Perez-Mercader, and 
M. Visser, Schwinger 's Dynamical Ca simir Effect: Bulk 
Energy Contribution, hep-th/9609195. Accepted for pub- 
lication in Physics Letters B. 

R. K. Pathria, Statistical Mechanics, (Pergamon, Oxford, 
England, 1972). 

R. Bahan and C. Bloch, Annals of Physics (NY) 60, 401- 
447 (1970). 

S. K. Blau, M. Visser, and A. Wipf, Nuclear Physics, 
B310 163-180 (1988). 
T. H. Boyer, Phys. Rev. 174 1764 (1968). 
^ , J T. H. Boyer, J. Math. Phys. 10 1729 (1969). 

yj(2^+ 1) [si{x)sii{x) + Sf^{x)s^{x)] — -x^, (AlO) [19] M. Abramowitz and I. A. Stegun, Handbook oi mathe- 

matical functions, (Dover, New York, 1972) 
I. S. Gradshteyn and I. M. Ryzhik, Table of integrals, 
series, and products, (Academic, New York, 1980) 



Note in particular that 



i=i 



[14] 

[15] 

[16] 

[17] 
[18] 
[19] 

[20] 



£(2^ + 1) W,ix)s'Ax) ~ s,{x)snx)] = -x', (All) 



and 



^(2£ + 1) [2seix)se{x) + s'eix)4{x) + se{x)s';{x)] 



= 2x\ 



(A12) 



20 



